Shell-model half-lives for the = 82 nuclei and their implications for the r-process 
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We have performed large-scale shell-model calculations of the half-lives and neutron-branching 
probabilties of the r-process waiting point nuclei at the magic neutron number A'' = 82. We hnd 
good agreement with the measured half-lives of ^^^Ag and ^"^"Cd. Our shell-model half-lives are 
noticeably shorter than those currently adopted in r-process simulations. Our calculation suggests 
that ^^°Cd is not produced in beta-flow equilibrium with the other N — 82 isotones on the r-process 
path. 

PACS numbers; 21.60.Cs, 21.60.Ka, 27.40. +z, 23.40.-s 



About half of the elements heavier than mass number 
A — 60 are made in the astrophysical r-process, a se- 
quence of neutron capture and beta decay processes [Q-^ . 
The r-process is associated with environment of rela- 
tively high temperatures (T 10^ K) and very high 
neutron densities (> 10^° neutrons/cm^) such that the 
intervals between neutron captures are generally much 
smaller than the P lifetimes, i.e. t„ <^ in the r-process. 
Thus, nuclei are quickly transmuted into neutron-richer 
isotopes, decreasing the neutron separation energy Sn- 
This series of successive neutron captures comes to a stop 
when the (n, 7) capture rate for an isotope equals the rate 
of the destructive (7, n) photodisintegration rate. Then 
the r-process has to wait for the most neutron-rich nuclei 
to /3-decay. Under the typical conditions expected for the 
r-process, the (n, 7) ^ (7, n) equilibrium is achieved at 
neutron separation energies, Sn ^ 2 MeV Q]. This con- 
dition mainly determines the r-process path, which is lo- 
cated about 15-20 units away from the valley of stability. 
The r-process path reaches the neutron shell closures at 
N — 50, 82, and 126 at such low Z-values that Sn is too 
small to allow the formation of still more neutron-rich 
isotopes; the isotopes then have to /3-decay. To over- 
come the shell gap at the magic neutron numbers and 
produce heavier nuclei, the material has to undergo a se- 
ries of alternating /3-decays and neutron captures before 
it reaches a nucleus close enough to stability to have Sn 
large enough to allow for the continuation of the sequence 
of neutron capture reactions. Noting that the /3-decay 
half-lives are relatively long at the magic neutron num- 
bers, the r-process network waits long enough at these 
neutron numbers to build up abundance peaks related to 
the mass numbers A « 80, 130, and 195. Furthermore the 
duration of the r-process, i.e. the minimal time required 
to transmute, at one site, seed nuclei into nuclei around 
A « 200, is dominantly given by the sum of the half-lives 
of the r-process nuclei at the three magic neutron num- 
bers. It appears as if the required minimal time is longer 
than the duration of the favorable r-process conditions 
in the neutrino-driven wind from type II supernovae , 



which is the currently most favored r-process site. 

Simulations of the r-process require a knowledge of nu- 
clear properties far from the valley of stability. As the rel- 
evant nuclei are not experimentally accessible, theoretical 
predictions for the relevant quantities (i.e. neutron sep- 
aration energies and half-lives) are needed. This Letter 
is concerned with the calculation of /9-decays of r-process 
nuclei at the magic neutron number N = 82. These (3- 
decays are determined by the weak low-energy tails of the 
Gamow- Teller strength distribution, mediated by the op- 
erator CTT— , and provide quite a challenge to theoretical 
modelling as they are not constrained by sumrules. Pre- 
vious estimates have been based on semi-empirical global 
models, quasiparticle random phase approximation, or 
very recently, the Hartree-Fock-Bogoliubov method. But 
the method of choice to calculate Gamow- Teller tran- 
sitions is the interacting nuclear shell model, and deci- 
sive progress in programming and hardware make now 
reliable shell model calculations of the half-lives of the 

= 82 r-process waiting point nuclei feasible. 

Our shell model calculations have been performed with 
the code ANTOINE developed by E. Caurier As model 
space we chose the 037/2, 1^3/2,5/21 2si/2, Oft.11/2 orbitals 
outside the = 50 core for neutrons, thus assuming 
a closed iV = 82 shell configuration in the parent nu- 
cleus. For protons our model space was spanned by 
the lpi/2, Ogg/2,7/2, 1^3/2,5/2, 2si/2 orbitals where a max- 
imum of 2 (3) protons were allowed to be excited from 
the lpi/2, O99/2 orbitals to the rest of the orbitals in the 
parent (daughter) nucleus. The Ogg/2 neutron orbit and 
the O/111/2 proton orbit has been excluded from our model 
space to remove spurious center of mass configurations. 
Therefore, we do not allow for Gamow- Teller transitions 
from the 0179/2 and ft.11/2 neutron orbitals, which should, 
however, not contribute significantly to the low-energy 
decays we are interested in here. The lpi/2 proton or- 
bit has been included to describe the l/2~ isomeric state 
seen in ^'^^In and expected in the others = 82 odd- A 
isotones, but does not play any role in the decay of the 
ground states. The residual interaction can be split into a 
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monopole part and a renormalized G-matrix component 
which can be derived from the nucleon-nucleon poten- 
tial. We use the interaction of ref. for the gdshii/2 
orbits and the KLS interaction |^ for the interaction be- 
tween the previous orbits with the orbit. To derive 
the appropriate monopole part, we followed the prescrip- 
tion given by Zuker joj, and fine-tuned the monopoles 
to reproduce known spectra of nuclei around the N = 82 
shell closure. As shell model studies overestimate the GT 
strength by a universal factor, we have scaled our results 
by the appropriate factor (0.74)^ pC[ |. 

The Qf3 values have been taken either from exper- 
iment (^^^In) or from the mass compilation of Duflo 



and Zuker Note that the Extended Thomas- Fermi 
with Strutinsky Integral approach (ETFSI) [|2| and the 
microscopic-macroscopic (FRDM) model of Moller 
give very similar Qp values (with typical uncertainties 
of 250 keV for 10 MeV) so that the associated 

uncertainty in the half-lives is small. 

The shell-model half-lives are summarized in Table || 
and are compared to other theoretical predictions in 
Fig. H For Z = 47-49, the half-lives are known 
experimentally and our shell model values are 

slightly faster. This, however, is expected as our still 
truncated model space will miss some correlations and 
hence slightly overestimates the Gamow- Teller matrix el- 
ements. 

Our shell model half-lives show significant and impor- 
tant differences to those calculated in the FRDM p^ ] 



and the ETFSI approach |12|, which have been typically 
used in r-process simulations. Although the latter pre- 
dicts a Z-dependence of the half-lives in the = 82 
isotones very similar to the present results, the ETFSI 
half-lives are longer on average by factors 4-5 indicat- 
ing that the method fails to shift enough Gamow- Teller 
(GT) strength to low energies |l2|. The FRDM half- lives 
show a very pronounced odd-even dependence which is 
predicted neither by ETFSI nor shell model. While the 
FRDM half-lives for odd-A A^ = 82 isotones approxi- 
mately agree with the shell model results (within a fac- 
tor of 2) and the experimental values for ^^^In and Ag, 
they overestimate the half-lives for even isotones by an 
order of magnitude. As such an odd-even dependence 
is not present in the experimental half-lives (nor in the 
r-process abundances) it is probably an artifact of the 
FRDM model. The absence of odd-even effects can be 
understood considering that the main contribution to the 
half- life comes from transitions from a gj /2 neutron to a 
gQ/2 proton due to the energy gap between the gg/2 and 
the other orbits. Therefore, neither the GT matrix ele- 
ments nor the half-lifes show a strong odd-even depen- 
dence along the A = 82 isotonic line. Noting that the 
values in the FRDM model are very similar to the ones 
used and that the main difference with our results ap- 
pears when the final nuclei is odd-odd we conclude that 
the odd-even effect must stem from the treatment of the 



pn interaction in the FRDM approach. Very recently, 
Engel et al. have performed half-life calculations of r- 
process nuclei within the HFB model . Unfortunately 
their studies are yet restricted to even-even nuclei only, 
but they obtain results which, except for a factor of 2, 
closely resemble the present shell model results. Ref. Q 
points out that the half-lives of the A^ = 82 waiting point 
nuclei are noticeably shorter than currently assumed in 
r-process simulations, in support of our findings. 

Odd-A nuclei in this mass range usually exhibit a low- 
lying 1/2^ isomeric state which can be related to a pro- 
ton hole in the Pi/2 orbital. These isomeric states can 
affect the r-process half-lives in two different ways: i) If 
low enough in energy, the isomeric state can be popu- 
lated thermally; ii) in a non-equilibrium picture the iso- 
meric state can be fed by the preceding neutron capture 
on the A^ = 81 nucleus. The half-live of the isomeric 
state has been measured in ^'^^In (350 ms), very similar 
to the ground state half- life (280 ms). We have calcu- 
lated the energy positions and half-lives of the isomeric 
states within our shell model approach. We find that the 
excitation energy of the isomeric state slowly decreases 
within the A^ = 82 isotones when moving from ^^'^Nb 
(£;* = 500 keV) to ^^^In (375 keV) where experimen- 
tally only the isomeric state in -'^^^In is known (at 360 
keV). Importantly our calculation predicts the half- lives 
of the isomeric states to be comparable to the ground 
state half-lives in all cases (see Table ||). Thus, the effec- 
tive r-process half-lives will be very close to the ground 
state half-lives. We note that the isomeric state in ^■^^In 
dominantly decays by first-forbidden transitions, as the 
approximately closed gg/2 proton configuration in this 
state strongly suppresses low-energy GT transitions. We 
calculate a half-live of the isomeric state of 274 ms due 
to first-forbidden decay, about 30% faster than the ex- 
perimental value (350 ms). For the A = 82 nuclei with 
Z < 47 the 59/2 orbital is not anymore closed for the 
isomeric state allowing for GT transitions at low-energy. 
Consequently these nuclei decay by GT transitions rather 
than first-forbidden ones. 

An interesting, but yet open question is whether the 
r-process proceeds in /?-flow equilibrium also at the wait- 
ing points related to magic neutron numbers If so, 
the duration of the r-process has to be larger than the 
sum of the beta half-lives of the nuclei in /?-flow equilib- 
rium. In this appealing picture, the observed r-process 
abundances scale like the respective /3-decay half-lives, if 
the former are corrected for /^-delayed neutron emissions 
during their decays from the r-process paths towards the 
stable nuclei which are observed as r-process abundances. 
Using our shell model /3-strength functions we have calcu- 
lated the probability Pi„ that the (3 decay is accompanied 
by the emission of a neutron, defined as the relative prob- 
ability of the /?-decay rate above the neutron emission 
threshold, Sn- For consistency we adopted the Sn values 
from Dufl-O and Zuker (DZ), which for a mild parity ef- 
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feet, gives similar results than the ETFSI model, while 
the FRDM model predicts a significantly slower decrease 
of Sn with decreasing Z. As the Pi„ values are rather 
sensitive to the neutron separation energies we have ac- 
counted for the differences between the ETFSI and the 
DZ predictions by assigning an uncertainty of 500 keV 
to the DZ Sn values. In this way we have calculated an 
equally probable range for the shell model Pi„ values as 
shown in Fig. |^. We find rather small neutron emission 
probabilities for Z = 46-49 (which will only slightly been 
increased if the 1/i-proton orbitals were included); for the 
smaller .Z- values our results approximately resembles the 
FRDM values, also showing a noticeable odd-even depen- 
dence. 

Using the solar r-process abundances |]l8| and the shell 
model neutron emission probabilities Pin, we have deter- 
mined the abundances of the = 82 progenitor nuclei, 
n(Z), on the r-process path by one-step iteration; we have 
checked that second-order branchings change the results 
only unsignificantly. If /3-flow equilibrium at the waiting 
points is indeed achieved, one has n{Z) ^ t{Z). Thus, 
up to a constant n{Z) can be expressed as the so-called 
/3-flow half- life Tfjf. Fixing the constant appropriately, 
figure ^ shows that /3-flow equilibrium can be attained for 
the = 82 isotones with Z = 44-47, but fails by more 
than a factor of 3 for ^^°Cd and ^'^^In. As the systematic 
of neutron separation energies put ^'^'^Cd on the r-process 
path, our results suggest that the conditions which allow 
to build the = 82 r-process abundance peak do not 
last long enough to achieve /3-flow equilibrium for this 
nucleus. This is consistent with the expectation that the 
r-process peaks at iV = 82 and N = 126 are made under 
different conditions Recent observations also 

indicate that the nuclides in the iV = 82 and 126 abun- 
dance peaks are produced at different sites . Further- 
more, the assumption of /3-flow equilibrium, which on the 
r-process path should be fulfilled the better the shorter 
the half- life, leads to unphysical (negative) T^/ values 
for Z < 43, indicating that these nuclei are not on the r- 
process path. However, firm conclusion here can only be 
reached after reducing the uncertainties in the Pi„ values 
for all nuclei in the decay sequence to stability. 

So far we have discussed the r-process as a sequence of 
competing neutron capture and /3-decay processes. If the 
r-process site is indeed the neutrino-driven wind above a 
newly born neutron star, then it occurs in a very strong 
neutrino fiux and charged-current {ve,^) reactions can 
substitute for (3 decays; in this case the picture of /3- 
fiow equilibrium has to be extended to 'wcak-fiow' equi- 
librium. McLaughlin and Fuller |2^ have shown that 
weak steady flow equilibrium at the iV = 82 waiting point 
normally cannot be attained in the neutrino-driven wind 
model. However, their study adopted /3 half-lives taken 
from jl^] which predict a Z-dependence of the half-lives 
in disagreement with the recent theoretical studies (in- 
cluding the present) and the data. Nevertheless, when 



we reinvestigate this question using the present shell- 
model /3-decay rates and the charge current neutrino 
rates of [|3|j2|], we flnd that ^^^Ag and "°Cd cannot be 
produce in /3-flow equilibrium within the neutrino-driven 
wind model in agreement with the conclusions reached in 
ref. H. 

In conclusion, we have calculated shell model half-lives 
and neutron emission probabilities for the N — 82 wait- 
ing point nuclei in the r-process, finding good agreement 
with the experimentally known half-lives for the Z — 47- 
49 nuclei. Our half-lives are significantly shorter than 
the ETFSI and FRDM half-lives, which are frequently 
used in r-process simulations. Our results indicate that 
i29^g jg produced in /3-fiow equilibrium together with the 
lighter isotones at the iV = 82 waiting point. R-process 
simulations usually include ^'^"Cd and even ^^^In in the 
r-process path at freeze-out. If so, they will not be syn- 
thesized in /3-flow equilibrium. That fact together with 
the shorter half-lives implies a shorter waiting time at 

= 82. This is quite welcome to remove possible con- 
flicts between the required duration time for the r-process 
and the expansion time scale of the neutrino-driven wind 
scenario which requires the r-process to occur in a frac- 
tion of a second. 
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FIG. 3. Test of /3-flow equilibrium by comparing the shell 
model half-lives for the TV = 82 r-process waiting-point nuclei 
with the /3-flow half-hves. The error bars in the latter reflect 
the errors in the Pi„ and in the solar r-process abundances. 
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as calculated in the FRDM and the present shell model ap- 
proach. 



TABLE I. Comparison of the shell model half- lives for 
the ground and isomeric state with experiment p"7p5| . All 
half-lives are in ms. 



Nucleus 


Ground state 


Isomeric state 




Expt. Theor. 


Expt. Theor. 


i^^In 


280 ± 30 177 


350 ± 50 274 




195 ± 35 146 






46 ± 9 35.1 


39.5 


128p^ 


27.3 




i^^Rh 


11.8 


12.4 


i^^Ru 


9.6 






4.3 


4.2 


Mo 


3.5 






1.8 


1.7 


1222, 


1.5 
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